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Abstract: We have realized a conical matter wave lens. The repulsive 
potential of a focused laser beam was used to launch a Bose-Einstein 
condensate into a radially expanding wavepacket whose perfect ring shape 
was ensured by energy conservation. In spite of significant interactions 
between atoms, the spatial and velocity widths of the ring along its radial 
dimension remained extremely narrow, as also confirmed by numerical 
simulations. Our results open the possibility for cylindrical atom optics 
without the perturbing effect of mean-field interactions. 
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1. Introduction 

Gaseous Bose-Einstein condensates (BECs) present unique possibilities as bright, coherent 
matter wave sources for applications to coherent atom optics, holography and interferome- 
try [1, 2, 3, 4, 5]. There is much interest in creating novel atom optical elements for focusing 
and deposition [6, 7], coupling of atoms with high efficiency into atom chip guides [8], col- 
limation of atomic trajectories within an atomic clock to suppress systematic effects [9], and 
for the exploration of fundamental physics with matter waves [10, 11]. However, unlike pho- 
tonic optics, propagation of matter waves is profoundly affected by the interactions between 
atoms. For example, in a gaseous Bose-Einstein condensate, the "mean-field" driven expan- 
sion destroys the Heisenberg-limited momentum distribution possessed by a condensate while 
it is trapped, creating a matter wave with a broad spread of velocities that is undesirable for 
atom-optical applications. Therefore, it is of great interest to discover methods of controlling 
and suppressing interaction effects during propagation of a matter wave field and within atom 
interferometers. Such approaches include manipulation of the scattering length using Feshbach 
resonances [12], Raman outcoupling of a BEC [13, 14], the use of extremely dilute condensates 
[15] and interferometry with fermionic gases where i-wave scattering is forbidden [16]. 

In this work, we report on the realization of a conical lens, or "axicon" for Bose-Einstein 
condensates (BECs) using the dipole force exerted by a focused, blue-detuned laser beam. To 
our knowledge this is the first observation of the far-field pattern of an axicon using neutral 
atoms. In our experiment, the atoms expand radially outward in the form of an extremely thin, 
circular ring in spite of significant interactions between atoms. This is in marked contrast to 
the conventional BEC expansion, which is characterized by strong spatial dispersion and a 
broad final velocity distribution. We have performed numerical solutions of the Gross-Pi taevskii 
equation confirming that this interplay between the external and internal forces can be exploited 
to tailor the expansion of a BEC using light fields. 

2. Experimental Description 

We produce Bose-Einstein condensates of roughly 3 x 10^ sodium atoms in a purely linear 
quadrupole potential formed by a pair of anti-Helmholtz coils following the procedure in our 
earlier work [17]. To create the BEC, we first use a blue-detuned "optical plug" consisting 
of a far-off resonance laser at 532 nm which has been focused to a 1 /e^ waist of 40 jJ.m to 
prevent Majorana transitions near the magnetic field zero [18, 19, 20]. The laser beam intensity 



is ramped to zero within 200 ms after the condensate is produced, resulting in a metastable 
condensate in a "linear" potential formed only by the quadrupole coils, with an axial field 
gradient of Z?t = 21 G/cm. The predicted radial and axial Thomas-Fermi radii of the condensate 
SiieRp =18 /im andR- —10 /im, respectively. 



Fig. 1. Artistic rendition of the operation of the atomic axicon. A condensate (wavefunction 
shown in black and white) is initially localized near the origin. The external potential is 
shown in false color, (a) At time r = 0, the linear potential is turned off, and a blue-detuned 
optical plug is turned on. (b) For times r > the atoms rapidly "roll down" the potential 
hill, forming a ring-shaped density distribution that expands radially outwards in the x — y 
plane. 

We created a ring-shaped matter wave output by simultaneously shutting off the quadrupole 
field and instantly turning on a second 532 nm laser beam focused to the atoms' location that 
rapidly expelled them. Hereafter we will refer to this expelling beam simply as "the plug beam". 
It propagated along the z-axis and had a 70 jxm radius (roughly 4 times larger than the 
cloud radius) and variable power The entire condensate was instantaneously located on the top 
of an azimuthally symmetric potential hill with a steep gradient and began to accelerate radially 
outwards, as illustrated schematically in Figure 1. This strong radial acceleration caused the 
atom cloud to form an extremely thin, nearly perfectly shaped ring that propagated outward 
(see Figure 2a for the ring structure and associated movie). 

To observe the expanding ring, we waited for a variable time delay after turning on the 
plug beam (t is hereafter denoted time-of-flight, or TOF). We then took an absorption image 
along the z-direction using a 200 microsecond probe pulse resonant with the transition from 
the F ^ I ground level to the F' = 2 excited level. After a few milliseconds, the atoms had 
been expelled from the region of the laser focus. After this point, neither the external potential 
nor the interactions played a significant role (the expansion had lowered the density to the point 
where interaction effects were negligible). The subsequent dynamics of the cloud was thereafter 
purely ballistic. 

We observed a narrow ring-shaped atom density distribution (see the movie in Figure 2a). 
Its width could only be measured accurately after carefully focusing the imaging system on the 
falling atom cloud (gravity is parallel to the z-axis). Since atoms gain exactly the same kinetic 
energy independent of which direction they were expelled, the nearly perfect circular shape was 
completely independent of minor misalignments of the plug laser or deviations from cylindrical 
symmetry of the laser focus. 

By contrast, the movie in Figure 2b shows the mean-field driven expansion of the BEC in the 
absence of the plug at a similar time-of-flight. The mean-field pressure creates a broad distribu- 
tion of radial velocities whose spread is of the order of the sound velocity in the condensate, 
about 1 cm/s. In the presence of the plug, however, this distribution becomes squeezed into a 
naiTow packet traveling at a large mean velocity of up to 6 cm/s. We observed such thin rings 
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Fig. 2. Far field pattern of a conical atom beam, (a) Movie of the ring expansion 
and (b) movie of ordinary BEC expansion [Multimedia available for download at 
http://www.physics.gatech.edu/chandra/index.htm]. The frame size is 2.6 mmx 2.6 mm. 
(c) Ring radius versus TOF and (d) velocity of expansion versus the peak AC Stark Shift, 
proportional to the laser intensity. The solid curve is a fit to a square root function arising 
from the energy conservation relation Eqn. 1 in the text. 



only for Bose-Einstein condensates, where the initial spatial and momentum distributions were 
simultaneously narrow. At our longest times-of-flight, these matter wavepackets had a spatial 
width which was 70 times smaller than the ring radius. Moreover, apart from an increase in the 
radius of the ring with time, the shape of the density distribution did not change substantially 
over the timescale of our observations. 

The ballistic nature of the expansion can be clearly seen in Figure 2c for times-of-flight 
greater than 5 ms. The maximum of the absorption occurs on a ring whose radius increases 
linearly with time. The slope of this line is simply the atomic velocity Vf , which we have 
plotted in Figure 2d as a function of the laser intensity, proportional to the peak AC Stark shift 
U, which was determined from the measured laser power and beam width. The velocity fits 
well to a square root function Vf °^ U'/^ with the expected prefactor for energy conservation. 

We can understand our observations by considering the focused laser beam as a refractive 
element for atom optics and by neglecting atom-atom interactions, at least initially. The inho- 
mogeneous dipole force shapes the transverse velocity distribution of the atoms in the following 
manner For a laser beam with intensity radius W the AC Stark shift potential is given by 
U{p) — Uoe^^P' . An initially stationary atom which is located at transverse position p sat- 
isfying < p < Rp, where Rp is the Thomas-Fermi radius of the condensate in the x — y plane. 



will be accelerated to an asymptotic radial velocity given by the energy conservation relation: 



Thus if Rp <C W, to lowest order each atom will be given the same positive radial velocity 



VfD = Y independent of its initial location. The correction terms will shape the final dis- 
tribution, as we will see later In an analogy with conventional optics, we identify the atoms' 
transverse velocity with the angle of propagation of light rays with respect to an optical axis. 
The coordinate along that axis is the equivalent of time in our experiment. In this picture, our 
atoms experience the equivalent of a negative axicon lens, which deflects incident light rays 
through a fixed angle independent of their transverse spatial position measured with respect to 
that axis. ' 

There have been relatively few reaUzations of refractive elements for atom optics. Most atom 
optics experiments have dealt primarily with thin lenses, where the change in trajectory during 
the interaction with the lens element can be neglected. Kaenders et al. have realized a posi- 
tive axicon for a slow atomic beam by exploiting the linear magnetic potential of a quadrupole 
magnet, which resulted in a uniform transverse force across the sample that caused the beam 
to focus on a line [21]. Since the laser is not pulsed in our experiment, but simply turned on 
instantaneously, the atoms' velocity distribution is significantly modified during the interaction 
with the beam. Thus we have realized a thick atom lens, which does not require the force to be 
constant in space as it was for the thin lens implemented by Kaenders et al. The conical feature 
of the final velocity distribution arises simply from energy conservation and an extremely nar- 
row spatial distribution of atoms Rp ^ W, and thus does not depend on the details of the forces 
which are applied. The implementation of thick lenses is greatly facilitated using Bose-Einstein 
condensates since the atomic dynamics are intrinsically slow. It is therefore straightforward to 
engineer a long interaction time with the lens element, which opens up new possibilities for 
realizing refractive elements for atom optics [22, 23]. 

The radius of the ring is determined by the peak AC Stark shift, which determines the en- 
ergy gained by atoms in traversing the plug potential. This is the equivalent of the cone angle 
of the axicon lens, which determines the deflection angle. However, the radial distribution of 
the atom cloud within the ring itself is more complex, and one must in general consider both 
the effects of interactions as well as the potential due to the lens. We have measured this dis- 
tribution for various times of flight and the data are shown in Figure 3 for our maximum plug 
Stark shift of kB x 6.5 jiK. Before application of the plug, the Thomas-Fermi density profile of 

the atoms is given by ni{p,z) = «o(l ^ —^z^), where «o is the peak density. After 

the interaction, the final density profile is given in cylindrical coordinates by «^(p, where 
(/) is the azimuthal angle. We measure the column density nf{p,(j)) ~ f nf{p ,z)dz- Due to 
the cylindrical symmetry of the problem, we can define an azimuthally averaged density pro- 
file nf{p) = 27tp / nf{p,<^)d<f), from which the normalization condition / nf{p)dp ~ N 



results, where is the number of atoms. 

In Figure 3a we have plotted the measured atom radial density for a succession of times- 
of-flight from 5 ms to 19 ms. The average allowed us to account for slight misalignments be- 
tween the plug and the Thomas-Fermi distribution of the atoms, which caused some azimuthal 

'it should be noted that there is one difference-the transverse velocity change is only constant for an initially 
stationary Bose-Einstein condensate, whereas an optical axicon is insensitive to the incident angle of the light rays as 
well as their position, with the change in angle being a constant in the paraxial approximation. For our atoms, this 
would con'espond to a non-zero radial momentum before the plug is applied. 
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variations in the atom density. Each 1 -dimensional data set was obtained by averaging the two- 
dimensional image using the following procedure. A 1 pixel-thin slice was taken through the 
image after it was rotated through an angle about the ring center, which could be readily de- 
termined to within the pixel size, about 5 jJ-m. For each slice, we folded the data about the center 
to obtain two data sets corresponding to angles 6 and 9 + n, both starting at p = 0. The pro- 
cess was repeated 180 times with 1 degree rotation increments and the final radial distribution 
for each image was the sum of all the data sets. For all times-of-flight, the distribution clearly 
has an asymmetric shape, with a long tail extending toward smaller radii and a relatively sharp 
cutoff at larger radii. We quantified the asymmetry by defining the half- width at half-maximum 
(HWHM) for the distribution with respect to its peak, defined as p = po. This allowed us to 
define a backward and forward HWHM, Whacn and Wf, respectively, corresponding to p < po 
and p > Po, respectively. These are plotted in Figure 3b, and clearly show that Whu^k > Wf for 
the majority of the data points. 
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Fig. 3. Spatial distribution of the ring. Data were taken at the maximum plug Stark shift of 
fcflX 6.5 /xK.(a) Radial number density of atoms for TOF from 5 to 19 ms in increments 
of 2 ms, showing the outward radial evolution, (b) Spatial Half-Width at half-maximum 
(HWHM) for forward width Wf (solid circle) and backward width Whack (square), (c) Time 
evolution of a classical "wavepacket" simulation of our experiment, as described in the 
text. The probability density is shown for times 0, 5, 10, 15 and 20 ms, and confirms the 
temporal focusing observed in (b). 



We can explain our results in terms of a simple atom-optical picture in which the plug laser 
acts as a lens that focuses the atomic distribution. To understand the physics, we concentrate 
on two groups of atoms-group A initially located at p close to and group B located at p 
close to Rp, the Thomas-Fermi radius of the BEC. At f = 0, atoms in group A experience a 
negligible force since the gradient in the laser intensity is close to zero. Therefore, they remain 
momentarily stationary. However, atoms in group B experience an immediate force due to the 



finite gradient, and begin moving radially outward. This leads to a momentary spreading of the 
wavepacket. This spreading is aided by atom-atom interactions, which act to reduce the atom 
density n further through the positive mean-field pressure °^ [24]. However, for longer times, 
group A atoms begin to move outward as well. They acquire more kinetic energy than group 
B atoms due to their larger initial potential energy, as expressed in Eqn. 1 . Eventually, group A 
will catch up with group B at some intermediate time Tf, which we denote the focusing time, 
when the spreading is reversed. This focusing time is, typically, considerably larger than the 
average transit time of atoms Xnans through the beam. The key is that for times Ttnms <t <Tf 
the atoms' velocity is close to the asymptotic value given by Eqn. 1, but their spatial distribution 
has not yet equilibrated. Moreover, by this time the axicon has reduced the atom density by 
spreading the atoms onto a ring, and interactions have no further effect on the propagation, 
which is governed purely by linear optics. Finally, at very long times t Tf, the wavepacket 
should defocus. That is, it will resemble a map of the asymptotic velocity distribution, with fast 
atoms in front, and slow atoms behind. 

Our experimental observations lie precisely in this intermediate timescale f ~ Tf. Thus we 
observe no increase in the forward width Wf and a reduction in the backward width Whack- 
The latter provides evidence of the temporal focusing of the wavepacket. Further evidence was 
provided by classical simulations which we have performed based on the specific parameters 
of our experiment. We numerically computed single atom trajectories in the presence 
of the plug for times < f < fi,„,. We chose various starting locations < Pi ^ < Rp at time 
t = 0, where k = 1,2, ...nmax, and computed a sequence of final locations Pf at f = fj,-,,,. The 
resulting classical "wavepacket" was constructed by assigning to each Pf a probability density 
weighted in proportion to the initial number of atoms located at p,jt, i.e.,o<: iTzpi f ni{pi i^,z)dz, 

the column density in the trap. The integration limits used are z ~ ±j^Rj, — p^, located on the 
Thomas-Fermi surface of the cloud. The result is shown in Figure 3c for = 200 trajectories 
and times fi,„, = 0,5, 10, 15,20 ms. The simulated wavepacket shows a strong resemblance to 
the data in Figure 3a. In particular, we observe an asymmetric distribution with a shrinking tail, 
i.e., Wbiick decreases with time. A close examination of the final classical wavepacket shows that 
it is in fact multi-valued at the leading edge, an indication that many classical trajectories have 
crossed one another, an effect confirmed by a detailed examination of the individual trajectories 
used to construct the wavepacket. At this point the classical analysis begins to break down, 
and one requires a full quantum treatment. Nonetheless, this simple simulation confirms the 
temporal focusing of the axicon with a time on the order of, or close to 20 ms. 

This "ray optics" picture cannot accurately be used to compute the atom distribution near 
the temporal focusing point. However, it is intuitively clear why the forward width Wf in our 
measurements remains stationary-atoms simply pile up from behind, and do not spread out 
until a much longer time t 3> Tf, when the different momentum components separate again. 
Unfortunately, we could not observe the ring at larger TOE as it exceeded the camera sensor 
size^. Although our experiments do not measure the asymptotic momentum distribution, we 
can make a rough estimate of the M^, or quality factor for our atom beam, which measures the 
phase space volume occupied. From the reduction of W/j^ck with time in Figure 3 we estimate a 
velocity width of order 5vr ^ 2 mm/s. The measured spatial width near the temporal focus is 
no larger than 20 jim, and therefore, = ^dr x dv,- < 30. While not yet at the diffraction 
limit, we note that typical values for RE outcoupled atom lasers are in the range of ~ 10, not 
far from our result [25]. 

In a further analogy with optics, our experiment may be described as a "cavity-dumped" 
atom laser, since the entire condensate is emitted as a single matter wave pulse. There are 

^Demagnification of the image would have allowed us to observe the ring at larger TOF, but at the cost of reduced 
spatial resolution, which would compromise the width measurement. 



close analogies between our work and a host of experiments that extract atoms slowly from 
a BEC, where the inhomogeneous mean field interaction between outcoupled atoms and the 
stationary condensate affects the spatial profile of the laser. For example, RF coupling results 
in a roughly constant energy but a broad transverse distribution due to the mean-field [26, 27], 
which can result in distortion of the spatial mode due to gravitational sag [25]. By contrast, our 
conical matter wave is extremely symmetric, and is relatively free of mode distortions due to 
the fact that gravity plays no important role in the extraction process. In this regard, our method 
is a significant step forward in comparison with RF techniques, especially when outcoupling 
away from the Thomas-Fermi surface. The conical mode shape could be useful for a number 
of atom-optical applications where the unique spatial mode can be exploited, for example, to 
form Bessel atomic beams that do not spread out in time. We plan to explore this avenue in the 
future. 





Fig. 4. Side view of expanding ring shown at 10 ms (top: a,b) and 15 (below: c,d) ms for 
both ring shaped expansions (left: a,c) as well as ordinary BEC expansion (right: b,d). Each 
image is 0.8 mm x 2.2 mm. 



While the plug created a narrow radial wavepacket, we found that the atoms' velocity along 
the z-direction was almost completely unaffected. Due to the long confocal parameter of 2.8 
mm, there is a negligible variation of the laser intensity and therefore no force exerted along 
that direction over the relevant experimental timescales. In Figure 4 we have observed the atoms 
after 10 and 15 ms TOF imaged along one of the radial directions. The ring of atoms appears 
in this view as a band with a considerable axial width, comparable to that of an ordinary BEC 
expansion in the absence of the plug, as shown in Figure 4(b,d). 

3. Theory 

While much insight was gained by a purely classical simulation, the full dynamics requires a 
quantum mechanical treatment. To this end, we have simulated the dynamics of the BEC in the 
presence of the plug in cylindrical coordinates p and z, taking advantage of the azimuthal sym- 
metry of the problem. We describe the condensate by the order parameter i//(p,z,f), normal- 
ized such that n{p,z,t) = A^|v/(p,z,f)p. The condensate dynamics is governed by the Gross- 
Pitaevskii equation 

m-^^xi/{p,z,t)^ (^^^V'- + V {p,z,t) +Ng\xi/{p,z,t)\"^ Y{p,z,t) (2) 

where is the number of atoms, g = AjtfP'a /m is the mean field interaction strength, m is the 
atomic mass, a is the scattering length, and V {p ,z,t) is the combined external potential applied 
by the magnetic trap and the plug beam. We assume the magnetic trap is turned off and the plug 
beam is turned on instantaneously at time f = 0. Explicitly, we write the potential as 



y(r,f) 



Uoexp 



2p^ 



,t<0 
,t>Q 



(3) 



where JJ-m ~ 27tti x 0.7 MHz/G is the atomic magnetic moment. Bp is the radial magnetic field 
gradient, Uo is the Stark-shift at the center of the plug, and W is the radius of the plug. At 
f = 0, we assume the condensate is in the ground state of the magnetic trap with the plug turned 
off. The ground state is itself calculated by evolving Eqn. 2 in imaginary time. The dynamics 
induced by the plug is then found for f > by integrating Eqn. 2 in real time. The grid size, 
spatial resolution, and time step were varied to ensure convergence. To ease the computational 
burden, the simulation was done for a smaller atom number than the experiment but with the 
same magnetic field gradient (see Table 1). We used approximately the same values of the ratio 
of plug potential to chemical potential and the ratio of plug radius to Thomas-Fermi radius 
of the condensate as in the experiment. 



Table 1. Comparison of theoretical and experimental parameters for ring shaped expansion. 
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Experiment 


70 


18 


6.4 nK 


0.7 11 K 


3 X 10^ 


11 



In Figure 5 we show the numerical results of the expansion, which provide key qualitative 
insights into our data. The initial density profile is shown in Figure 5a, and closely follows the 
Thomas-Fermi distribution. After applying the plug potential for a time 1.3 ms, the simulation 
was stopped and the final density distribution computed as shown in Figure 5b. Figure 5c shows 
the corresponding radial momentum distribution of the wavefunction. By this time the density 
had been lowered by a factor of 200 due to the expansion caused by mean-field repulsion of 
atoms as well as the radial force of the plug potential. Moreover, at this point in time, the peak 
of the atom density has moved beyond the 1/e^ radius of the plug potential. Therefore, the 
subsequent expansion of the cloud should be mostly ballistic, with only small corrections to the 
momentum distribution expected. 

The density profile in Figure 5b shows that the cloud has expanded considerably in the z- 
direction. Moreover, in the radial direction, the atoms are being forced outward and form a 
steep front at p = 32 jim. However, a spatial tail of the distribution still exists for smaller 
p. Our data in Figure 3a for 5 ms and at short TOF bears a striking similarity to the spatial 
distribution in Figure 5b. Figure 5c shows the velocity distribution of the atoms and confirms 
that the majority of the atoms occupy a very narrow peak in momentum space near 38 mm/s, 
corresponding to a kinetic energy of 1.95 jlK, about 90% of the peak plug potential. This is 
consistent with our observations in Figure 3b of an extremely narrow velocity spread for the 
expanding ring. The velocity width (HWHM) computed from the simulation is about 600 jim/s, 
about a factor of 3 larger than the Heisenberg limited width of 215 jim/s for a cloud of radius 
Rp = 6.4/xm. For comparison we have also shown the momentum distribution of a mean-field 
driven expansion in the absence of the plug. To understand the role that mean-field effects play 
on the ring formation, we performed similar calculations for an expansion time of 1.3 ms using 
the same initial wavefunction but with the mean-field term in the Gross-Pitaevskii equation set 
to zero. This represents the purely "linear" evolution of the atomic wavepacket. The resulting 
momentum distribution is also shown in Figure 5c, and is qualitatively very similar to the non- 
linear evolution. A similar sharp peak with a narrow velocity width is formed, demonstrating 



that the dynamics are largely insensitive to the effects of interactions. 

Two other features emerge from the simulations-in both the interacting and non-interacting 
cases in Figure 5c there is a long tail extending to low velocities, most likely due to the slow 
evolution of atoms initially located very close to p = 0. These atoms will catch up with the peak 
of the wavepacket only at later times. Moreover, it is intriguing to note the presence of smaller, 
"satellite" momentum peaks in both simulations which appear to be interference phenomena. 
To explore these further, we extended the non-interacting calculation to the longer time of 1.8 
ms. This calculation proved easier than the fully interacting case due to the larger grid size 
used. Due to the similarity between the curves in Figure 5c, the linear evolution provides a 
useful guide to the dynamics. An expanded view of the momentum peaks is shown in Figure 
5d at both 1.3 and 1.8 ms for the linear evolution. The interference peaks appear to increase in 
amplitude with time, and coincide with a reduction of the low velocity tail. It is clear that the 
radial phase profile is not flat, and there may be differences between the interacting and non- 
interacting cases at longer times. Future work will explore these interferences in detail, with 
an eye toward understanding their effect on atom-optical applications of the axicon lens, where 
the phase profile is of importance. 



Fig. 5. Numerical simulation of ring shaped expansion. Computed density profiles in the p- 
z plane for (a) the initial Thomas-Fermi distribution and (b) at 1.3 ms, after the atoms have 
mostly traversed the region containing the plug. In (c) we show the velocity distribution of 
the wavefunction of (b) (solid black) as well as that obtained for a purely linear temporal 
evolution, i.e., by turning off the nonlinear term in the Gross-Pitaevskii equation during 
the expansion (dot-dashed red). For comparison, we also show the velocity distribution 
obtained by mean-field expansion without the plug (dashed blue). In (d) we examine the 
region near the velocity peak for the purely linear expansion at times 1.3 (dot-dashed red) 
and 1.8 ms (solid blue). For each curve, the velocity distribution is normalized to its peak 
value. 
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4. Conclusion 



In summary, we have demonstrated a conical lens for atoms that produces an extremely narrow 
velocity distribution, which could have applications to interaction-free matter wave manipula- 
tion. To our knowledge, this is the first observation of the far field pattern of such a conical lens 
in atom optics. Our results point to the power of using Bose-Einstein condensates, where one 
can readily apply thick refractive elements to shape the spatial and momentum distributions. 
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